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摘要 ”纳米 通道 内 的 反 电 渗析 过 程 是 一 种 将 纳米 通道 两 端 由 于 浓度 差 造 成 的 能 量 势 差 直 接 转 化 成 电能 的 新 方法 。 由 于 反 
渗析 中 局 部 浓度 随 通道 位 置 变化 ， 壁 面 电荷 密度 也 会 随 通道 位 置 发 生变 化 。 本 文 引 入 双 电 层 带电 Basic-Stern 模型 ， 将 
壁面 电荷 密度 与 局 部 浓度 关联 ， 并 在 传统 空间 电荷 模型 的 基础 上 ， 预 测 了 反 电 渗析 过 程 的 电流 -电压 曲线 和 电动 势 。 修 正 
的 空间 电荷 模型 由 于 引入 双 电 层 带电 模型 ， 将 难以 测量 的 物理 量 ， 如 壁面 电荷 密度 ， 转 化 为 较 易 测 得 的 物理 量 和 已 知 的 
模型 参数 。 与 过 去 模型 的 对 比 表 明 修 正 的 空间 电荷 模型 在 合理 选择 模型 参数 的 条 件 下 更 符合 实验 结果 。 
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Abstract: Reverse electrodialysis is a kind of energy conversion method that generates electricity from the free 


energy in salinity gradient. Due to concentration gradient along the streamwise direction, surface charge density 


will consequently vary along channel. In this work, we introduce the Basic-Stern model to relate the surface 


charge density with local concentration and predict current-voltage curve and electromotive force in the reverse 


electrodialysis, based on the conventional space charge model. Our modified model converts the surface charge 


density, which is difficult to measure, into accessible physical quantities and model parameters. Compared with 


the previous models, our modified model fits the experimental results better with proper parameters. 
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0 引 言 
盐 差 能 是 不 同 浓度 的 溶液 自发 混合 时 可 以 利用 由 于 吸附 或 分 解 过 程 带电 ， 从 而 在 溶液 中 形成 双 电 
的 自由 能 差 ， 广泛 存 在 于 江河 入 海口 处 。 反 电 渗 析 层 。 双 电 层 内 异种 电荷 离子 〈counter-ion ) 浓度 远大 
是 一 种 将 盐 差 能 转换 成 电能 的 方式 。 利 用 离子 选择 于 同 种 电荷 离子 〈co-ion) 浓度 。 由 于 纳米 通道 的 尺 
作 膜 的 传统 反 电 渗析 过 程 在 1954 年 首次 由 Pattle"] ”十 与 双 电 层 厚度 相当 ， 甚 至 远 小 于 双 电 层 厚度 ， 异 
提出 并 实现 。1976 年 ，Weinstein 和 Leitz 中 设计 了 一 种 电荷 离子 将 在 通道 内 占 主 导 。 当 通道 两 端 存在 浓 
种 电池 组 装置 应 用 于 反 电 渗析 ， 并 且 给 出 了 优化 操 度 差 时 ， 由 于 双 电 层 效应 ， 纳 米 通道 呈现 出 对 不 同 
作 条 件 的 理论 依据 。 近 年 来 ， 随 着 膜 技 术 的 发 展 ， 电 性 离子 的 选择 性 ， 从 而 产生 净 电 和 荷 的 传输 。 构 成 
反 电 渗析 重新 得 到 重视 ， 大 规模 的 应 用 也 进入 试验 回路 后 ， A 
阶段 申 。 除 此 之 外 ， 由 于 纳米 技术 的 快速 进步 ,无 机 图 1 为 Kim 等 四 设计 的 纳米 通道 反 电 渗 析 系 统 
材料 制 成 的 纳米 通道 被 应 用 于 反 电 渗析 中 ， 为 小 尺 的 示意 图 ， 采 用 Ag/AgCl 作为 极 板 。 实 验 表 明 ， 纳 
度 下 能 量 转化 提供 了 一 种 新 的 可 能 四。 ed a 
当 电 解 质 溶液 与 固体 接触 时 ， 固 体 表面 自发 的 出 1-3 个 量 级 。 同 时 ,无 机 纳米 通道 稳定 性 更 好 ,不 
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会 像 有 机 膜 一 样 缩水 和 膨胀 。 除 此 之 外 ， 依 赖 于 成 
熟 的 CMOS 流程 的 纳米 通道 加 工 技术 相 较 于 传统 的 
膜 工 艺 成 本 也 较 低 外 。 最 近 ，Kim 等 站 已 研究 出 纳米 
通道 阵列 的 反 电 渗析 系统 ， 大 大 提高 了 电流 密度 ， 
纳米 通道 阵列 组 装 形成 的 无 机 膜 将 有 可 能 大 大 拓展 
反 电 渗析 的 应 用 范围 。 
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1 纳米 通道 反 电 渗析 示意 图 


Fig.1 Schematic for reverse electrodialysis in nanochannel 


已经 有 了 初步 成 果 ， 但 是 对 于 纳米 通 

道内 反 电 渗析 过 程 的 理论 研究 仍然 较 少 。Fair 和 
Osterlelg 最 早 采用 空间 电荷 模型 (space-charge mod- 
el) 研究 了 反 电 渗析 过 程 ， 但 是 其 研究 针对 的 是 有 机 
膜 内 的 圆柱 形 纳米 通道 ， 并 未 应 用 于 无 机 纳米 通道 
的 分 析 ， 另 外 ， 模 型 中 采用 了 给 定 壁 面 电荷 的 假设 。 
Kim 等 由 采用 简化 的 基于 Teorel-Meyer-Sievers 
(TMS) 模型 的 一 维 模型 分 析 了 电动 势 和 迁移 数 等 
参数 ， 得 到 了 定性 上 正确 的 结果 ， 但 是 在 通道 尺度 
很 小 或 者 溶液 池内 浓度 较 高 时 存在 较 大 误差 。Kim 
等 中 和 Guo 等 加 分 别 采 用 有 限 差 分 和 有 限 元 方法 计 
算 了 整个 系统 的 电场 和 浓度 场 ， 他 们 均 采 用 给 定 壁 


尽管 实验 


由 电 蓓 密度 的 边界 条 件 ， 没 有 考虑 浓度 变化 对 电荷 
密度 的 影响 。 之 前 的 研究 已 经 表明 [5]， 表 面 电 荷 密 
度 依 赖 于 局 部 浓度 等 因素 。 因 此 ， 给 定 壁 面 电荷 密 
度 的 假设 在 纳米 通道 反 电 渗 析 系 统 中 难以 成 立 。 男 
外 ， 当 双 电 层 厚 度 大 于 通道 尺度 时 ， 双 电 层 重 膨 导 
致 的 离子 富 集 效应 3 中 也 不 可 忽略 。 


1 模型 与 方法 


考虑 两 个 无 穷 大 平板 组 成 的 二 维 纳米 通道 ， 平 
板 间 距 h 即 为 通道 高 度 。 在 TMS 模型 中 ， 忽 略 电势 
和 浓度 在 截面 方 铝 的 变化 ， 仅 考虑 轴 问 的 
Nernst-Planck 方程 : 


Wy =W(X) (2) 


以 及 进出 口 处 的 Donnan 平衡 。 这 里 ，c 为 i 离子 的 


c= (x) 
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浓度 ， Ce D, ，z; 分 别 为 i 离子 的 扩散 系数 
和 电荷 数 ， @、 了 分别 为 元 电荷 量 ，Boltzmann 
常数 和 温度 。 es T 取 室温 298 K。 
在 空间 电荷 模型 中 ， 将 电势 作 如 下 分 解 : 
D=G (x) +y (x,y) (3) 


根据 Gouy-Chapman 的 双 电 层 模型 ,y 满足 截面 
上 的 Poisson-Boltzmann (PB) 方程 


2 2 Y 
2 A | 4 | (4) 
€ kT 
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其 中 ，c (x) 是 截面 上 的 表征 浓度 ， 当 双 电 层 不 重印 
时 ，c (x) 即 为 通道 中 间 浓 度 ， 当 双 电 层 重 县 时 ， 

c (x) 视 为 外 推 至 无 穷 远 处 的 浓度 。 根 据 Basic-Stem 
模型 5 可 以 得 到 不 同 的 c (x) 时 的 壁面 电荷 密度 : 

kT —O kTInl0 o 

“(c) In (pH—pK) 


(5) 
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在 给 定 o 的 边界 条 件 下 ， 通 过 求解 不 同 截面 处 
的 PB 方程 ， 得 到 不 同 截 面 处 的 电势 分 布 。 在 偏离 平 
衡 态 不 太 远 的 条 件 下 , 考虑 如 下 关于 电流 工 和 溶质 通 
量 J 的 线性 唯 象 方程 组 多 ， 


dlnc” dy 

I=k,| -RT [| Rr | (7) 
dlnc” dy” 

J =k,| -RT 了 js kz 坚 | (9) 


其 中 ，k,i,j=1,2 是 唯 象 系数 。 

通过 求解 PB 方程 得 到 电势 分 布 ， 在 低 Re 数 的 
条 件 下 ,流动 方向 速度 可 以 由 不 可 压 的 Navier-Stokes 
方程 得 到 ，u =u(y,c (x),f "(x)) 。 通 过 积分 截面 不 
同位 置 的 通 量 , 可 推出 电流 了 的 表达 式 , 与 唯 
组 比较 后 即 得 到 k ,i,j=1,2 的 形式 ,在 二 维 纳米 通 


内 ， 以 有 为 例 : 


1 ae se 
LnInaxIV 呈 1 


ele 合 。 计算 结 果 对 pH, pK 值 比较 敏感 ， 调 整 pH, pK 值 

Ej; 壬 | cosh¥ (H? -yjdy 可 以 得 到 比较 接近 实验 的 结果 ， 但 是 难以 保证 反 电 
7 0 

2wc pH ( 即 曲 线 的 横 截 距 和 和 斜率 )。 这 可 能 是 进出 口 的 
1 (Dp(-Y)+ DexpYY) 和 y 9) Donnan 平衡 条 件 不 能 计 全 满 中 导 至 的。 实际 上 ， 对 
sw(c 本 于 高 度 26 nm 的 通道 ， 按 照 进口 高 浓度 和 出 口 低 浓 

+ J, coshy| fo sinh Ydydydy 度 计算 的 德 拜 长 度 分 别 为 9.2 nm 和 29.2 nm, 进口 双 
电 层 未 完全 重 闯 ， 而 出 口 己 完 全 重 疙 ， 这 可 能 对 结 


参考 Sasidhar 等 ("的 数值 方法 ， 首 先 将 式 (8) 除 。 ” 果 带 来 影响 。 
以 式 (7), 得 到 J 的 表达 式 , 从 通道 入 口 至 出 口 积 4 


0.016 - 


求解 积分 方程 得 到 J 。 然 后 将 是 (8) 代 入 式 (7)， 消 il es 
SS A ee 
去 dg 》 积分 后 可 求 得 了 。 eT O、 ~ -spKe75 HES.6 
dx 0010F ~ 以 、 -< 人 -pkK=7.7 pH=5.6 
全 oo 京 、 ->-pK=7.9 pH=5.6 
疆 \, y = 0.008 上 上 ~ Ro 和 ~ 一 太一 Experiment 
3 一 口 果 与 讨论 EF te ~ 尽 、 
呈 . 、 ee 3 ~ 
Gn 模拟 中 采用 Kim 等 "的 实验 参数 ， 电 解 质 溶液 a a 
为 KCl,， 浓 溶液 端 浓度 Cr=103 molL， 稀 溶液 端 浓 0002[ Rg 
度 C=10“ mol/L， 纳 米 通道 长 宽 高 分 别 为 140 jm， oop | , I 
25 hm 和 26 nm。 如 图 2 所 示 , TMS 模型 的 模拟 结果 8 
与 实验 存在 较 大 的 偏差 ， 同时， 由 于 TMS 模型 中 的 Ee 


壁面 电荷 密度 09 是 未 知 量 ， 只 能 采用 拟 合 的 万 式 与 。 图 3 不 同 pH, pK 时 修正 空间 电荷 模型 得 到 的 电流 -电压 曲线 
实验 对 比 。 


Fig.3 Current-voltage curve from modified space charge 
model for different pH, pK 
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图 2 不 同 电荷 密度 下 TMS 模型 得 到 的 电流 -电压 曲线 40| ， 本 
107 107 10 
Fig.2 Current-voltage curve from TMS model for different sur- G, molD 
G=10” (mol/l) 


face charge densities 


采用 修正 的 空间 电荷 模型 计算 的 结果 如 图 3 所 全 加 

示 ， Basic-Stern 模型 中 的 参数 如 下 方式 选取 : C og 002C/m, 人 的 空 | ne MR 
尿 实 验 中 当 Im 六 

厂 =8 nm?; pH 采用 Kim 等 的 实验 值 中 ， 由 于 实验 误 ee 
差 ，pH =5.6+0.3 ， 由 于 实验 并 未 给 出 不 同 条 件 下 浓度 比 ， 而 改变 通道 高 度 时 ， 修 正 的 空间 电荷 模型 
的 pH 5 本 文中 采用 实验 误差 范围 内 人 允 F 的 值 与 实验 得 到 的 结果 比 TMS 模型 和 必 系 始 空 | 间 电荷 模 型 更 接近 
进行 拟 合 ， 根据 Hiemstra 等 的 表面 杂 化 模型 “， 取 ”实验 值 。 图 4 和 图 5 中 取 pH=5.3，pK=7.9 ， 在 前 述 
反应 常数 PK = 7.5 , 但 是 考虑 到 本 文中 可 能 出 现 双 电 ”的 误差 范围 内 与 实验 结果 最 接近 
层 重 个， 同时 pK 作为 界面 的 本 征 参数 ， 其 具体 大 小 
仍 存在 争议 ， 本 文中 取 pK =7.5~ 7.9 与 实验 进行 拟 
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中 


TMS 模型 、 空 间 电 和 荷 模型 、 修 正 空间 电荷 模型 得 
到 的 不 同 浓度 比 时 的 电动 势 


Fig.4 Electromotive force from TMS model space charge model, 


modified space charge model for different concentration ratios 
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5 TMS 模型 、 空 间 电 荷 模型 和 修正 的 空间 电荷 模型 得 
到 的 不 同 通道 高 度 下 电动 势 


Fig.S Electromotive force from TMS model space charge model, 


modified space charge model for different channel heights 


4 结 论 


日 于 修正 模型 综合 考虑 了 壁面 带电 的 复杂 物理 
机 制 ， 建 立 了 壁面 电荷 密度 与 局 部 浓度 的 关联 ， 可 
以 通过 易 测 的 物理 量 精确 计算 得 到 难以 测量 得 到 
壁面 电荷 密度 ， 而 不 需要 将 壁面 电荷 密度 作为 拟 合 
参数 与 实验 比拟 。 同 时 ， 计 算 结 果 表 明 ， 修 正 的 空 
间 电 荷 模型 在 取 合适 参数 的 条 件 下 ， 比 TMS 模型 和 
未 修正 的 空间 电荷 模型 与 实验 结果 符合 得 更 好 。 
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